Subdiffusive model of substance releasing from a thick membrane 
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We study both theoretically and experimentally the process of subdiffusive substance releasing 
from a thick membrane. The theoretical model uses the subdiffusion equation with fractional time 
derivative and specific boundary conditions at the membrane surfaces. Using a special ansatz we 
find analytical formulas describing the time evolution of concentration profiles and an amount of the 
substance remains in the membrane. Fitting the theoretical functions to the experimental results, 
we estimate the subdiffusion coefficient of polyethylene glycol 2000 in agarose hydrogel. 

PACS numbers: 05.40.-a,66.10.-x 
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Subdiffusion qualitatively differs from the normal dif- 
fusion. It occurs in a medium where mobility of parti- 
cles is strongly hindered due to internal structure of the 
medium, as for example in porous media or gels [11, Q ■ 
The subdiffusion is characterized by the relation where 
the mean square displacement of a Brownian particle is 
a power function of time 2] 



2D^ 



r(i + 7) 

is the subdiffusion coefficient measured in the units 



the theoretical functions with our experimental results 
we estimate the subdiffusion coefficient of polyethylene 
glycol 2000 {PEG2Qm) in agarose hydrogel. 

The system consists of three homogeneous parts which 
are separated from each other by two infinitely thin 
partially permeable membranes located at x = xi and 
X = X2 (see Fig. [T]). In each part there are the same 
subdiffusion parameter 7 and the subdiffusion coefficient 
(1) Dj. In the following these parts will be denoted as 1 



m^/s^ and 7 is a subdiffusion parameter which obeys 
< 7 < 1. The case of 7 = 1 corresponds to the normal 
diffusion. Till now, there have been only a few methods 
to extract the subdiffusion parameters from experimen- 
tal data (see for example [l|, Q)- We mention here the 
method of measuring the time evolution of near mem- 
brane layers in a system with one thin membrane [l[. 
The method utilized the model of the subdiffusive trans- 
port in the system where a thin membrane separates pure 
solvent from homogeneous solution. 

In this paper we present the model of releasing of the 
substance from a thick membrane. The system under 
consideration is assumed to be homogeneous in a plane 
perpendicular to the x axis which is perpendicular to the 
membrane surfaces. Thus, the system is effectively one- 
dimensional. The thick membrane is treated here as a 
homogeneous slab limited by two thin membranes. We 
find the theoretical functions describing a concentration 
of the transported substance and the time evolution of 
the amount of the substance which remains in the thick 
membrane. This model can be used to extract the subd- 
iffusion parameters from experimental data. Comparing 



for X < xi, M for Xi < x < X2 and 2 for a; > a;2. We 




FIG. 1: The schematic view of the system under considera- 
tion. The concentration in the interval (a, b) is equal to the 
initial one Co. 

consider the system where its middle part is filled with 
a homogeneous solution and the external parts contain a 
pure solvent at the initial moment. 

Let us assume that the transport process is described 
by the subdiffusion equation with the Reimmann- 
Liouville fractional time derivative 



•Electronic address: 'tkoszt@pu.kieIce.pl I 
t Electronic address: dworeckiQpu. kielce.]^ 
■1 Electronic address: |kale@amg.gd.pl| 



dC{x,t) 



d^C{x,t) 



(2) 



dt ''dt^-r dx^ ' 

where C{x,t) denotes the concentration of transported 
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substance. The initial condition is 

{0 , X < xi , 
Co , xi <x<X2 , (3) 
, X > X2 ■ 

To solve the subdifFusion equation in the three-part 
system one needs six boundary conditions. Two of them 
demand vanishing of the solutions at ±cxd 

Ci(-oo,t) = 0, C2(oo,t)=0, (4) 

two others demand the continuity of the fluxes at the 
membrane surfaces 

Jl{x^ ,t) ^ Jm{x^ ,t) , Jm{x2 ,t) = J2{X2 ,t) , (5) 

where the subdiffusion flux is given by the formula 
J,{x,t) = -D^d^-f /dt'^-''dQ{x,t)/dx, i = 1,M,2. The 
main assumption of our model is that the missing bound- 
ary conditions at the thin membranes are fixed as 

Ci{x^,t) = \i{t)CM{xt.t) , 

C2{xt,t) = \2{t)CM{x2,t) , (6) 

where 

Ai(0 =ai-6ie-"'^* , = aa - 626""'^*, (7) 

oi, a2, 61, 62, wi and W2 are positive constants. The 
functions ^ were found by trail and error to fit the ex- 
perimental results. The parameters \i{t) and \2{t) con- 
troll the permeabilityof the membranes; their interpreta- 
tion is presented in [j| . We note that previously we used 
the boundary conditions (H])-® where Ai and A2 are in- 
dependent of time ([4,] and references cited therein). 

It is easy to see that the exact analytical solutions 
of Eq. with the initial condition ^ and the bound- 
ary ones dH)-® are extremely hard to obtain when Ai 



and A2 are given by ([7]). To facilitate the problem we 
adopt the following assumption: when t <C tg, where 
tg is the average time when a particle passes the dis- 
tance between the thin membranes, the solutions of the 
system in near membrane regions can be obtained as 
for the system with one membrane. The motivation of 
this assumption is that for sufficiently small time a par- 
ticle localized in the vicinity of a membrane 'does not 
feel' the presence of another membrane. According to 
the above assumption, the concentration in the middle 
part of the system changes only in the relatively small 
near-membrane intervals, (xi,a) and (6,2:2) shown in 
Fig. ([T]). The points a and b play only auxiliary and 
illustrative role in our considerations and we do not con- 
sider their exact localization which actually changes in 
time. The concentration for x £ (a, b) remains unchanged 
Cm{x, t) — Cq. The presence of the interval (a, b) in the 
experimental plots suggests that the above assumption 
makes sense. The parameter tg can be estimated from 
the relation ^ putting d ^ X2 ~ xi = (Ax'^it)) which 
gives tg = (d2r(l-7)/2Z?^)i/7. 

Despite the fact that above assumption makes the cal- 
culations simpler, the solutions of Eq. ^ are still hard to 
obtain. To find approximate solutions we use the follow- 
ing ansatz: we solve Eq. (0) with the initial condition ^ 
and the boundary ones ^-(S^ for Ai and X2, which are 
independent on time, and next we replace Xi and A2 with 
Xi{t) and X2{t) (Eq. respectively, in the obtained 

solutions. 

Taking into account the above assumptions, we get 
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1 + Ai(t) ■'-1,7/2 



, X e {xi,a) 



CM{x,t) = Co , 



Co 



A2(t) f 

1 + A2(t) •'-1,7/2 



X G (a, b) 



, X e (6,2:2) 



(9) 



where 

^,p(^;«)-^|: fc,r(J,_,) {-^y ■ (11) 

(\ The function /i/,p(t;a) can be expressed in terms of the 

^ 2^ — 2:2 j ^^Q^ Pox .//-function 5]. It is easy to check that the functions 

J ' ©-(dni) obey the conditions 0-® and they approxi- 
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FIG. 2: The scheme of the experimental setup, see the text 
for a more detailed description. 



matcly solve the subdiffusion equation We mention 
here that the numerical solutions performed with the pa- 
rameters extracted from the experimental data and given 
in the further part of this work are very close to the func- 
tions ([8|l- p^ (this problem will be discussed in details 
elsewhere Jij). 

The experimentally measured function, which is more 
frequently used then the concentration profiles, is the 
time evolution of the amount of substance released 
from the sample R(t) 7]. In our system R{t) = 
J^^^Ci{x,t)dx + C2{x,t)dx; so, the amount of the 
substance, which remains in the sample i?M(i), is equal 
to RM{t) = Cod — R{t). After simple calculations from 
Eqs. (Ill) and (fTO| we get for in the long time approxima- 
tion 



RAiit) - Co 



d- 



1 + Xi{t) 



l + X2{t)J r(l-f7/2) 



(12) 



We apply our theoretical model to describe the releas- 
ing process of P_EG2000 from the agarose hydrogel. The 
measurement has been conducted in a membrane sys- 
tem shown in Fig. [21 The membrane system under study 
is a cell with three glass cuvettes separated by horizon- 
tally located membranes. Initially, we fill the lower and 
upper cuvettes with the agarose hydrogel solvent while 
in the middle cuvette there is an aqueous gel solution 
of transported substance. Then, the substance diffuses 
from the middle cuvette to the exterior ones through the 
membranes. Since the concentration gradients are in the 
vertical direction only, the diffusion is expected to be one- 
dimensional (along the axis x). The substance concen- 
tration is measured by means of the laser interferometric 
method ^ . The experimental set-up was already used 
to study transport in a system with one thin membrane 
and it is described in detail in the papers @. Here we 
only mention that it consists of the cuvette with two thin 
membranes, the Mach-Zehnder interferometer including 
the He — Ne laser, TV — CCD camera, and the comput- 
erized data acquisition system. For each measurement 



we prepared two gel samples: the pure gel 2% (w/v) wa- 
ter solution of agarose and the same gel dripped by the 
solute of PEG2000. The concentration of solutes in the 
gel was fixed to be 0.0075 mol/dm?. The agarose gel 
water solvent was prepared by dissolving agarose powder 
(Sigma) in 90'^C water. All experiments were performed 
at room temperature (22 ± 0.5)"C. The agarose gels are 
assumed to be inert to the solute at our experimental 
conditions. The polymer membranes (which are of the 
thickness 20 /j,to) initially separated the homogenous gel 
solution in one cuvette from the pure gels in another 
ones. At the beginning of the experiment the cuvettes 
were pressed to each other in close contact so that the 
diffusion across the membranes was initiated. For tech- 
nical reason the measurement of the concnentrations can 
be perfomed in the one part of the system only. 

In Fig. [3] we present the experimentally measured con- 
centrations in the middle part of the system for several 
times from to 7200 s. The errors of the concentra- 
tions are estimated as 10% of its value. The subdiffusion 
parameter 7 = 0.86 ± 0.03 was found in another exper- 
iment when the time evolution of the near membrane 
layer was analyzed by means of the method presented in 
. The theoretical functions, which also show in Fig.[3l 
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FIG. 3: The concentration profiles for the times given in the 
legend. Symbols represent the experimental data, continuous 
lines represent the theoretical functions. For clarity of the plot 
the error bars and the concentrations in the interval (2.5, 7.5) 
are not shown (inside this interval the experimental curves 
contained the constant function C = Co are present for all 
times) . 

are calculated for Co — 0.0075 mol/dm^, 3.1 x 

10^1" Ai(t) = 0.307 - 0.292 exp(-t/3043.5) 

and X2{t) = 0.309 - 0.287 exp (-t/2804.5). The subdif- 
fusion coefficient D^, which is independent of time, was 
treated as a fit parameter which ensures the best match- 
ing of theoretical and experimental results. The func- 
tions Ai (t) and A2 (t) were founded in the following man- 
ner. For each time the values of Ai and A2 give the best 
fit of the function © to the experimental data. Next, 



4 




1000 2000 3000 4000 5000 6000 7000 8000 
f 



FIG. 4: Time dependence of Ai and A2 obtained as fit pa- 
rameters (symbols) and their approximation given by Eq. ([7]) 
(continuous lines) for the parameters presented in the text. 



the values of Ai_2, presented in Fig.|4l were fitted by the 
functions ([7]) using the least square method. 



In Fig. [5] we present the time evolution of the amount of 
the substance which remains in the thick membrane. The 
experimental values were calculated for the data given 
in Fig. [3] by means of the numerical integration and the 
theoretical ones were obtained from Eq. (fT2|) . We observe 
a very good agreement of theoretical and experimental 
functions. 



Our analysis allows one to extract the subdiffusion co- 
efficient of the releasing substance from the experimental 
data. For PEG2000 transported in 2% agarose hydrogel 
we find = (3.1 ± 0.9) x 10"^° m^/s"-^^. Let us note 
that for t <^ l/cji ^ as well as for t ^ l/a;i,2 the func- 
tions Ai and A2 can be approximated by the constant 
functions, then R{t) ^ f^^. Thus, measuring the time 
evolution of the substance released from the membrane 
one can find the subdiffusion parameter 7. 



The parameters extracted from experimental data give 



to 



10^. It confirms that the functions 



TUl) can be 
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used to model the concentrations for times used in the 
experiment. When t is of the order of tg (or larger) we 
expect that Ai and A2 are constant, then the analytical 
solutions with the constant ratio of substance concentra- 
tions on both sides of the membrane surface can be used 
to describe this process (see (H). 



FIG. 5: The amount of the substance which remains in the 
thick membrane. Symbols represent the data obtained from 
the experiment, the continuous line represent the function 
(|12l) for the parameters given in the text. The error bars, 
which are estimated as 10% of values Rm, are not shown in 
the plot. 
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